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A cosmology calculator that computes times and distances as a func- 
tion of redshift for user-defined cosmological parameters is available at 
http://www.astro.ucla.edu/~wright/CosmoCalc.html. This note gives the for- 
mulae used by the cosmology calculator and discusses some of its implementa- 
tion. A version of the calculator that allows one to specify the equation of state 
. parameter w and w' is at ACC.html, and a version for converting the light travel 
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£j '. 1. INTRODUCTION 



times usually given in the popular press into redshifts is at DlttCalc.html. 
Subject headings: cosmology: miscellaneous - methods: miscellaneous 



There are many definitions of distances in cosmology, and it is often frustrating to 
calculate the kind of distance needed for a given situation. In particular, press releases 



(Clavin, 2006) very often use the light travel time distance, Dm = c(t — t em ) } even though 
this distance has very many undesirable characteristics. I wrote the cosmology calculator in 
order to be able to calculate Dm along with the observable angular and luminosity distances 
Da and Dl- It also computes the proper radial distance which is the only distance that is 
compatible with the common definition of distance for moving objects: the spatial separation 
at a common time. The common time used is now, or t . This it is called D now below. D now , 
also known as the comoving radial distance, is not measurable from a single position, but it 
can in principle be measured by a family of comoving observers. It is the one distance that 
satisfies the Hubble Law exactly: dD/dt = H(t)D without any corrections. Other distances 
are defined using quantities taken from more than one time and generally only satisfy the 
Hubble Law to first order. 
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2. EQUATIONS 

The following results have been given many times before. Hogg (1999) gives equations 
equivalent to the ones in this paper. But for completeness in documenting the cosmology 
calculator the equations are spelled out below. 

The metric is given by 

ds 2 = c 2 dt 2 - a(t) 2 R 2 (d^ 2 + S 2 (^)[d6 2 + sin 2 6d(t> 2 }) (1) 

where S(x) is sinh(x) if Q tot < 1 and sin(x) for Q tot > 1. The radius of curvature is given 
by R = (c/H )/ a/ 1 1 — Vltot | since I use the normalization that a = 1 at the current time so 
the redshift is given by 1 + z — 1/a. The past light cone is given by cdt = a{t)R dip so the 
comoving radial distance is 

„ „ , f cdt f 1 cda . . 

D now = R ij= —= — (2) 

J a Ji/(i+z) aa 

and of course the light travel time distance is given by 

D ltt = [cdt= f — = c{t - t{z)) (3) 

J Jl/{l+z) a 

The exact dynamics of the Universe can be captured using the the energy equation from the 
Newtonian approximation to cosmology, which gives 

a 2 GM(< a) 

= const. (4) 

2 a w 

GR modifies the acceleration equation by including the pressure as a source of gravity, but 
this just cancels the variation of the enclosed mass M(< a) caused by the pressure, so the 
energy equation is exact. So we can write a as H a \^X with 

X(a) = tt m /a + tt r /a 2 + Vt v a 2 + (1 - Q tot ) w 1 + 2q z + ... (5) 

where q is the deceleration parameter. Let us define a quantity Z that is the comoving 
distance divided by the Hubble radius cj H Q . Then 

f 1 da f z dz 9/ . . 

Z= - 7 == / — «*-* 2 (l + go )/2 + ... (6) 

J 1/(1+2) ay X Jo (l + z)VX 

To first order this agrees with the redshift, so Z = z + . . . This definition gives D now = 
(cZ/H ) and Dm = (c/ H ) J 1 /( 1+2 ) da/y/X. Then the angular size distance, which is defined 
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as D A — R/9 where R is the transverse proper length of an object that subtends an angle 9 
on the sky, is given by: 

c S(y/\l-Q tot \Z) 



D A = 



H (l + z)^\l^i 



tot\ 



= -^^(i + -(i-n tot )z 2 + —(i-n tot ) 2 z 4 + ..) (7) 

(1 + *) V 6 12(T ; J K J 

We can define a function J{x) given by 

sin(v / — x)j \J— x, x < 0; 

J(x) = <( sinh(Vi)/ v /x, x > 0; (8) 

1 + x/6 + x 2 /120 + ... + x n /(2n + 1)! + . . . , x « 0. 



Then 



da 

#o 7i/(i+ 2 ) a/X 
c f 1 da cZ 
Ho Ji/(i+z) a\fX H Q 



n n J([l^ojg) 

1 + Z 

D L = (l + z) 2 D A (9) 

Note that the luminosity distance -D L is defined to make the inverse square law work perfectly 
for bolometric fluxes, so that F bol = L/(A-nD 2 L ) for an object of luminosity L. 

The enclosed volume calculation requires the integral of either sin 2 i\)di\) or sinh 2 tpdip. 
The closed case becomes 

V-(< 2 ).4^I^^(|_^ (10) 
Note that ip = — VL to t\Z(z). The open case gives 

n<;) = 4^f(f^-f) (ii) 

The ratio of V(< z) to the Euclidean formula (47r/3)(i? ^) 3 is given by 



V (< Z ) n V- 6 X 4 " U, O \ *7 { \ 



21 n 



(4^/3)^)3 - tt(2n + 3)! 



1 + ^[(1 " n«)S(*) 2 ] + ^[(1 - n tot )Z(z) 2 ] 2 + ... (12) 



in both the open and closed cases. 
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3. IMPLEMENTATION 

The cosmology calculator is implemented as a Web page (CosmoCalc.html) that has 
a large number of Javascript definitions in the header, followed by immediately executed 
Javascript that writes a frameset to the current document. The frameset calls for the input 
form (CCform.html) and the output page (CCout.html). If Javascript is not enabled, or if 
there is an error in the Javascript, then the body of CosmoCalc.html is displayed. This body 
is just an error message saying that Javascript must be enabled. 

I have received several requests for the code used in the cosmology calculator. But 
since the code is in Javascript, it is in the HTML files in ASCII form. It is easy to save 
the page source to get the code, and it is easy to modify the code using any text editor. 
If your modifications introduce an error, then you will see the error message saying that 
Javascript must be enabled. This only means that you must find the error in your modified 
CosmoCalc.html. I have had to do this dozens of times, so don't be discouraged. 

Even if you do not intend to modify the code, downloading the three HTML files will 
let you run the calculator locally when not connected to the Internet. 

The numerical evaluation of the integrals for Z and t is done using the mid-point rule 
with 1000 panels. While this is not a very efficient use of CPU time, it was very easy to 
code. And with Javascript, the CPU time is consumed on the user's computer instead of 
the server. The functions being integrated go smoothly to zero as a goes to zero as long as 

n r > o. 

Another hidden aspect of the cosmology calculator is that it automatically sets the 
radiation density Vl r h 2 to the value appropriate for T D = 2.72528 K and three massless 
neutrino species, Q r h 2 = 4.165 x 10~ 5 . Here h = H o /100 km/sec/Mpc, and this factor 
includes a small (< 1%) boost in the neutrino density due to a slight transfer to e + e~ 
annihilation energy into neutrinos (Hannestad & Madsen, 1995). So if you want to verify 
some simple cases, like the empty Universe, you should use a large value of the Hubble 
constant which reduces the relative importance of this radiation term. For example, the 
open button with H Q = 97.78 and Vt M = gives an age of the Universe of 9.934 Gyr which 
is 0.7% from the expected 10 Gyr. But using H Q = 977.8 gives an age of 999 Myr which is 
only 0.1% from the expected 1 Gyr. 

For very early times, the age of the Universe is adjusted to allow for the existence of 
e + , e~, fj, + ,fi~, etc. using the g* and g*s for the standard model of particle physics given by 
Kolb & Turner (1990). 
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4. Light travel time inversion 



A slightly modified version of the cosmology calculator is at 
http://www.astro.ucla.edu/ ~wright/DlttCalc.html. 

The input form asks for the light travel time in Gyr instead of redshift. The redshift is found 
by evaluating the integral for light travel time in steps of —0.001 in a starting from a = 1 
until the input value is exceeded, and then interpolating to get a and thus z. Once z is 
known the calculations proceed as in CosmoCalc.html. 



Since the cosmology calculator was first written, there have been two developments that 
change the kinds of models that people want to run. I have created a new version of the 
cosmology calculator with more options and placed it in 
http : // www. astro, ucla. edu / ~ wright / ACC . html 

One development is the discovery of neutrino oscillations, indicating that the assumption of 
massless neutrinos is not correct. The neutrino temperature is (4/ll) 1 / 3 T Q = 1.95 K, and a 
typical momentum for a thermal neutrino is ~ 3/cT/c. This corresponds to the rest energy 
of a neutrino with mc 2 ~ 0.0005 eV. Since the neutrinos thermalized while still relativistic, 
their distribution is that of a relativistic Fermi-Dirac particle, so the neutrino density is given 



The number of spin states for a neutrino is one, but to allow for anti-neutrinos one should 
use g s = 2. 

The integral over neutrino momentum has to be evaluated for each step of the inte- 
gration over a so it needs to be done efficiently, even when the work is done on the user's 
computer. In the low temperature limit when kT u « m u c 2 the integral over x evaluates to 
m„(3/4)£(3)r(3), while in the high temperature limit kT v » m u c 2 the integral evaluates to 
(/cTj,/c 2 )(7/8)C(4)r(4). Both limits can be evaluated correctly by approximating the integral 
using 
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MORE OPTIONS 



by 




(13) 



^/m 2 u + x 2 (kT v ) 2 /c A 



x 



2 dx 



(3/4)C(3)r(3) v 7 ! + {x x kT v lm v c 2 ) 2 (14) 



exp(a;) + 1 



m 



•v 
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with the single integration knot at x 1 = [(7/8)C(4)r(4)]/[(3/4)C(3)r(3)] = 3.151. This 
approximation has a maximum error of < 3%. But a better approximation is 

/ y/ml + x*(kT v y/c* f d * - = m,(3/4)C(3)r(3)/( a ; 1 A;T,/m i ,c 2 ) (15) 
J expi^xj + i 

with 

/(y)«(l + /) V/3 (16) 
and P = 1.842 which has a maximum error < 0.3%. The mass at which neutrinos are 
semi-relativistic at the current epoch is 

m re ic 2 = x x hT v = 0.000531(T o /2. 72528 K) eV (17) 

The final result is that the effective neutrino density can be written as 

T \ 3 Ee^^/KeKl + Z)fm v ) 



^ h = U^SkJ ' 93.14 eV/e* (18) 

The normalization of 93.14 eV is from Mangano et al. (2005) and is 1.05% higher than the 
nominal due to residual coupling of annihilation energy to the neutrinos. In the relativistic 
limit the density is 1.53% higher. Increasing T u by 0.48% over the nominal (4/ll) 1 / 3 T when 
computing m re i allows for the relativistic limit, and this boost is included in Eq(17). 

This Q v (z) is just a function of z which gives Q v at z = and should not be confused 
with the actual Vl u at z's other than zero. In ACC.html the neutrino density at z = is 
subtracted from the input Om giving separate Qcm = — Q u (z = 0) for the CDM and 
baryons, and fl u for the neutrinos. Of course the neutrinos are not included in the radiation 
so a r h 2 = 2.778 x 10" 5 (T o /2.72528) 4 . 

For a hierarchical neutrino mass pattern, with masses 0.001, 0.009 & 0.049 eV, the 
change in distances introduced by neutrino masses is negligible. At redshifts up to 5 in the 
WMAP concordance flat ACDM cosmology with Vl M = 0.27 and H — 71, the changes are 
less than 0.01%. This mass pattern is the default when ACC.html is loaded. But even the 
more massive nearly degenerate neutrino mass patterns such as ~ 0.33, 0.33 & 0.33 eV have 
a minimal effect on the distances and times. 

The second change in the paradigm is the introduction of the equation of state w(z) = 
P/pc 2 as a parameter in the model. I have implemented w(z) as w = w + 2w'(l — a) 
following Linder (2003) who added the factor of "2" normalization to the simple form used 
by Chevallier & Polarski (2001). This functional form behaves well at high redshift and it 
allows an analytic solution for the dark energy density as a function of z. This solution is 

Pde = Pde,o a ^ 3+3wM ^ exp(-6u/[l - a]) (19) 
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Fig. 1. — The 4 distance measures Da, At*, D now and D L from bottom to top plotted 
against a fifth distance measure: the redshift z, for 3 different cosmological models: an open 
CDM model, a flat ACDM model, and a closed CDM model. The center and right panels 
are consistent with the CMB angular power spectrum, while the left and center panels are 
consistent with large scale structure. Only the center panel is consistent with supernova 
data. 

The defaults are w — — 1 and w' — and in that case ACC.html reduces to the Cosmo- 
Calc.html case. 

Unlike the neutrino masses, changes in w have substantial effects on distances and ages. 
Changing w to —0.7 instead of the default — f changes the age of the Universe by 6% and 
the luminosity distance at z = 1 by 7% when VL M = 0.27, Vt DE = 0.73, and H = 71 are 
left unchanged. But for w = —0.7 the model consistent with both supernovae (Riess et al. 
2004 gold+silver) and the CMB (Page et al. 2003 peak positions) changes to an open model 
with Vl M = 0.19, VL DE = 0.74 and H Q = 84, and this gives an age change of 17% and a 
luminosity distance change of 21% at z = 1. However, if Qm and Qde are allowed to vary 
as free parameters, the observable supernova signal in D L is reduced to only 0.3% (Wright 
2005). 

Finally ACC.html allows one to input T Q , allowing for easier tests of simple cases. 
With these changes, the X function is 
X(a) = (n CM + £l v (z))/a + Q r /a 2 + Q DE a -(i+3t»<>+6«,') ex p(-6u/[l - a]) + (1 - Q tot ) (20) 
The rest of the calculation of distances is unchanged. 
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6. CONCLUSION 

These cosmology calculators are suitable for interactive use, providing fairly quick an- 
swers for single cases. Users who wish to use the code for large scale calculations should 
translate it into a compiled language and change the quadrature formula. In particular, if 
fitting to datasets with redshifts and distance, the data should be sorted by redshift and 
the distance integrals evaluated for all objects in a single pass through the sorted data. 
Plotting figures is an obvious case where redshifts are computed in order: Figure 1 shows 
the distance measures discussed in this paper as a function of z for 3 different models. The 
formulae presented here were translated into Postscript for this figure, resulting in a 4 kB 
file. 
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